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With the Cowling approximation, the torsional oscillations on relativistic stars in tensor-vector- 
scalar (TeVeS) theory are examined. The spectrum leatures in TeVeS are very similar to those in 
general relativity (GR), but the torsional frequencies in TeVeS become larger than those expected 
in GR. We find that, compared with the fluid oscillations with polar parity, the torsional frequencies 
depend strongly on the gravitational theory. Since the dependences of fundamental frequencies on 
the gravitational theory and on the equation of state are different from those of overtone, it could 
be possible to distinguish TeVeS from GR in the strong-field regime via observations of this type of 
■ oscillations with the help of the observation of stellar mass. 
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I. INTRODUCTION 



It is well-known that the theory of general relativity (GR) is valid in a weak gravitational field such as our Solar 
System, as has been shown via many experiments. Looking at the case in the strong gravitational field, the situation 
is quite different from the case in the weak- field regime, i.e., the tests of gravitational theory in the strong- field 
regime are still very poor. However, with the development of technology, it is becoming possible to observe compact 
objects with high accuracy, and with these observations it will be possible to test the gravitational theory in a strong 
gravitational field [l|. In practice, one could observe compact objects not only via X-rays and 7-rays but also via 
> , gravitational waves emitted from the objects. So far, some possibilities for distinguishing the gravitational theory 
in the strong-field regime have been suggested. For instance, it has been proposed that one could distinguish the 
i—i, scalar-tensor theory [2| from GR by using the surface atomic line redshifts [3] or gravitational waves radiated from 
the neutron stars [J] . The possibility of a definitive test for GR with the direct observation of gravitational waves has 
r ^ \ also been pointed out @. 

■ As an alternative gravitational theory tensor-vector-scalar (TeVeS) theory has attracted considerable attention in 
recent years, as proposed by Bekenstein @ to extend the modified Newtonian dynamics @,|8[ into a relativistic theory, 
^\ , i.e., a covariant theory. As well as the modified Newtonian dynamics, TeVeS can explain the galaxy rotational curve 
and Tully-Fisher law without the presence of dark matter [6]. This theory is also successful in explaining strong 
, gravitational lensing Q and galaxy distributions through the evolving Universe without cold dark matter [10j |. For 
the strong gravitational region of TeVeS, Giannios found the Schwarzschild solution [ll|, Sagi and Bekenstein found 
the Reissner- Nordstrom solution [T^, and Lasky et al. produced the static, spherically symmetric stellar models in 
I ■ Te VeS by deriving the Tolman-Oppenheimer- VolkofF (TOV) equations in TeVeS [l3| . Additionally, Lasky and Donova 
examined the stability and quasinormal oscillations of the Schwarzshild solution in TeVeS [14j . 

Furthermore, some have suggested how to distinguish TeVeS from GR observationally. For instance, one can reveal 
the gravitational theory in a strong gravitational field with the redshift of the atomic spectral line radiating from 
the neutron star surface (l3j . with Shapiro delays of gravitational waves and photons or neutrinos [l5| . with the 
spectrum of gravitational waves emitted from compact objects 16], and with the rotational effect of neutron stars 
17]. In fact, via observations of gravitational waves due to stellar oscillations, one can obtain information about the 
stellar parameters, such as mass, radius, rotational rate, magnetic fields, and equation of state (e.g., [l8l - l22l ]). which 
is called "gravitational wave asteroseismology." In addition to this, the direct detection of gravitational waves could 
possibly be used to determine the radius of the accretion disk around a supermassive black hole [23[ or to know the 
magnetic effect during a stellar collapse [24[ . The test of gravitational theory in the strong- field regime is also one of 
the significant benefits of directly detecting gravitational waves. 

In this article, we focus on the torsional oscillations in neutron stars, which are incompressible oscillations. Usually, 
the frequency of this type of oscillations should degenerate to zero if the neutron stars consist of matter without 
elasticity. However, realistic stellar models have the solid crust region near the stellar surface, and in this region shear 
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torsional oscillations could exist. The boundary between the solid crust region and the fluid core is still uncertain 
because that depends strongly on the nuclear symmetric energy [25[ , but the density of this boundary is proposed to 
be 2.4 x 10 14 g/cm 3 [2f| or 1.28 x 10 14 g/cm 3 '27]. In this article, we adopt the density of 2.0 x 10 14 g/cm 3 to see the 
dependence of torsional frequencies on the gravitational theory. From the observational point of view, the observed 
quasiperiodic oscillations in giant flares [28j are considered to come from this type of oscillations [29j. 

In order to see the dependence of torsional frequencies on the gravitational theory, we adopt the relativistic Cowling 
approximation as a first step in this article, i.e., the perturbations except for the fluid one are neglected. It is known 
that this approximation is quite good at least in GR, because the torsional oscillations do not involve the density 
variation. A more detailed study including the perturbations of other fields will be done in the near future. 

This article is organized as follows. In the next section, we mention the stellar model in TeVeS, and in Sec. IIIII we 
derive the perturbation equations for torsional oscillations with the Cowling approximation. In Sec. IIV1 the concrete 
calculations of the obtained perturbation equations will be performed to show the dependence on the gravitational 
theory. Finally, we conclude in Sec. |Vj In this article, we adopt the unit of c = G = 1, where c and G denote the 
speed of light and the gravitational constant, respectively, and the metric signature is (— , +, +, +). 



II. STELLAR MODELS IN TEVES 

In this section, we only mention the fundamental parts of the theory. (See @ for details of TeVeS.) TeVeS is based 
on three dynamical gravitational fields: an Einstein metric g^u, a timelike four- vector field U^, and a scalar field <p, 
in addition to a nondynamical scalar field a. The vector field fulfills the normalization condition with the Einstein 
metric as g^vli^W = — 1, and the physical metric g^ is defined as 

9^ = e' 2ip g^ - 2U^U V sinh(2(^). (2.1) 

All quantities in the physical frame arc denoted with a tilde, and any quantity without a tilde is in the Einstein 
frame. The total action of TeVeS, S, contains contributions from the three dynamical fields mentioned above as well 
as the matter contribution jfjj, which includes two positive dimensionless parameters, k and K, corresponding to the 
coupling parameters for the scalar and vector fields. The field equations for the tensor, vector, and scalar fields can 
be obtained by varying the total action with respect to g^ v , , and <p, respectively. (See @ for the explicit forms.) 
Since a previous study for the neutron star structures in TeVeS has shown that the stellar properties are almost 
independent from the scalar coupling k [13], in this article we focus only on the dependence of vector coupling K. 
The restriction on K has not been discussed in great detail in the literature, but the authors of[l3| mentioned that 
K has to be less than 2 to construct the stellar models, and also that K should be less than 1 to produce a realistic 
stellar mass. According to this restriction on K, in this article we examine K as varying in the range of < K < 1. 
It should be noticed that GR is corresponding to the special case that K = and k = in TeVeS. 

As a background stellar model, we consider nonrotating relativistic stars, which have been investigated in [l3j]. A 
static, spherically symmetric stellar model can be expressed with the following metric: 

ds = -e v dt 2 + e c dr 2 + r 2 (d9 2 + sin 2 6d<j) 2 ) , (2.2) 

where v and C are functions of the radial coordinate r. Although the vector field on a static, spherically symmetric 
spacetime can be generally described as U.^ = (U t ,U r , 0, 0), for simplicity we adopt in this article the simple case that 
the radial component should be zero, i.e., W = 0, which is the same assumption adopted in [H, [H, G3- Then, with 
the normalization condition, the vector field can be fully determined as = [er v ^ 2 1 0, 0, 0). With this vector field, 
the physical metric is 

ds 2 = -e» +2 *dt 2 + e^dr 2 + r 2 e - 2 ^ (d9 2 + sin 2 Bd<t> 2 ) , (2.3) 

and the fluid four- velocity is u M = e^U^. Regarding the stellar matter, we assume a perfect fluid described by the 
energy-momentum tensor 

f^v = (p+ p) u^Uv + pg^ v , (2.4) 

where p and p are the pressure and energy density in the physical frame, respectively. Furthermore, in order to 
construct the stellar model, one needs to prepare the relation between p and p, i.e., the equation of state (EOS). In 
this article, we adopt the same EOS's as in j4|], which are polytropic ones derived by fitting functions to tabulated 
data of realistic EOS's known as EOS A and EOS II. With these EOS's, the maximum masses of a neutron star in 
GR are M = 1.65M for EOS A and M = 1.95M for EOS II. That is, EOS A and EOS II are considered as soft and 
intermediate EOS's, respectively. At last, the stellar modes in TeVeS can be constructed by using the recipe shown 
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III. TORSIONAL OSCILLATIONS 



As mentioned above, in this article, we focus on the torsional oscillations with the relativistic Cowling approximation, 
where we consider only the fluid perturbation with axial parity and the other perturbations of vector and tensor fields 
are neglected. It should be noted that the perturbation of scalar field exists only for polar perturbation. Considering 
the torsional oscillations, the Lagrangian displacement vector for the fluid perturbation can be expressed as 
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(3.1) 



where dg denotes the partial derivative with respect to 6, while Pi = Pi(cos9) is the Legendre polynomial of order I. 
Then, the nonzero component of the perturbed fluid four- velocity in the physical frame can be written as 



5u* = e-^^dtZ^- 



■dePt. 



(3.2) 



where d t denotes the partial derivative with respect to t. The perturbed energy-momentum tensor including the 
contribution from shear is given by 



STfj, v = (p + p) {du^Uu + u^riu) - 2fi8S^ u , 



(3.3) 



where fx is the shear modulus and is the shear tensor defined via = CnS^ 30, 31]. Here a^ v is the rate of 
shear tensor, which is defined as 



where P^ is the projection tensor 
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The speed of shear waves can be expressed as ^ = fi/(p + p) [3l[ , where a typical value of v s is around 10 8 cm/s in the 
crust of neutron stars. In this article, we adopt the shear modulus inside the crust determined by this simple relation. 
Finally, one can get the equations describing the fluid perturbations by taking a variation of the energy-momentum 
conservation law, i.e., i5(V pT a ^) = 0, which reduces to \7p5T al3 = with the Cowling approximation. The explicit 
form with a — d> becomes 
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where ( ) and (') denote the partial derivative with respect to t and r, respectively. Assuming that the perturbed 
variable has a harmonic time dependence, such as Z(t,r) = Z(r)e lult , Eq. (|3.6[) reduces to 
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Imposing appropriate boundary conditions on this equation, the problem to solve becomes the eigenvalue problem. 
In practice, the above equation will be integrated only in the crust region, i.e., the boundary conditions are imposed 
at the basis of crust (r = R c ) and at the stellar surface (r = R), because fi = in the fluid core as mentioned in the 
Introduction. In this article, we impose a zero traction condition at r = R c and the zero-torque condition at r = R 
[3ll |. Both conditions correspond to Z' = 0. 



IV. OSCILLATION SPECTRA 



In this section, we examine the torsional oscillations in the crust region of neutron stars both in GR and in TeVeS. 
With respect to the torsional modes in GR, it is known that one can see the dependence of £ in the fundamental 
modes, while the frequencies of overtone are almost independent of £, as shown in the left panel of Fig. [TJ where 
the specific frequencies in GR are plotted as a function of the stellar mass. On the other hand, we have done the 
numerical calculations to determine the frequencies in TeVeS with different value of K . As an example, the result 
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FIG. 1: Frequencies of fundamental modes and first overtone are plotted as functions of the ADM mass with EOS A in GR 
(left panel) and in TeVeS with K = 0.5 (right panel). In the figures, the solid lines correspond to the fundamental frequencies 
with £ = 2, 3, 4, and 5, i.e., (£, n) = (2, 0), (3,0), (4,0), and (5,0), where n is the number of nodes in the eigenfunctions, while 
the broken lines correspond to the first overtone with I = 2, 3, 4, and 5. 



with K — 0.5 is shown in the right panel of Fig. [JJ From this figure, it is shown that the feature of torsional modes 
in GR can be kept even in TeVeS. 

In order to see the dependence of frequencies of torsional oscillations on the gravitational theory, we focus on the 
1 = 2 oscillation modes. Figure [5] shows the fundamental frequencies (left panel) and the frequencies of first overtone 
(right panel), in both GR (solid lines) and TeVeS (broken lines) for the stellar models with EOS A. Similar to Fig. [TJ 
the frequencies are shown as functions of the stellar mass. From this figure, one can easily observe that the frequencies 
expected in TeVeS are quite different from those in GR. In practice, depending on the value of K, the frequencies in 
TeVeS become 34% larger for the fundamental oscillations and 150% larger for the first overtone than those in GR. 
In Fig. [3j we draw a figure similar to Fig. [2] but EOS II. One can observe that the relative augmentation of the 
frequencies in TeVeS compared with those in GR becomes 28% for the fundamental oscillations and 116% for the 
first overtone. That is, the stellar models with softer EOS seem to be more sensitive about the gravitational theory 
than those with stiffer EOS. Additionally, it should be emphasized that, compared with the fluid oscillations of polar 
parity accompanied by variation of density, the frequencies of torsional modes depend strongly on the gravitational 
theory, where the fre quen cies of polar parity in TeVeS become around 20% larger than those expected in GR 
And, as mentioned in |l6| . since this deviation between the frequencies in GR and in TeVeS results from the existence 
of a scalar field, observing the torsional oscillations as well as the fluid oscillations of polar parity could tell us the 
existence of the scalar field. 




FIG. 2: Frequencies of fundamental modes and first overtone with 1 = 2 are plotted as functions of the ADM mass with EOS 
A. Left and right panels correspond to the fundamental and first overtone frequencies, respectively. In both figures, the solid 
and broken lines correspond to the frequencies in GR and TeVeS, respectively. 

Furthermore, in Fig. |4j we plot the fundamental frequencies (left panel) and the first overtone (right panel) of 
torsional modes as functions of parameter K, where the ADM masses are fixed to be 1AM®. Here, the frequencies 
with I = 2, 3, 4, and 5 are shown in the left panel, while those with only £ = 2 are shown in the right panel, because 
the frequencies of overtone are almost independent of £, as we saw in Fig. [TJ Additionally, in both panels, the results 
for EOS A and EOS II are plotted using the solid lines with circle marks and the dotted lines with square marks, 
respectively, and the results in GR are shown at K = 0. From these figures, one can observe that the qualitative 
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FIG. 3: Similar to Fig. [2] but with EOS II. 



dependencies of frequency on the value of K are independent of the adopted EOS and the value of I. That is, the 
frequencies of torsional modes are increasing as the value of K becomes large. From the left panel of Fig. since 
the dependence of K on the fundamental frequencies is almost comparable to that of EOS, it might be difficult to 
distinguish the gravitational theory by using only the observation of fundamental modes. On the other hand, the 
right panel of Fig. 2] shows that the frequencies of overtone depend strongly on EOS rather than on the value of K . 
So, with observation of the frequencies of overtone, it could be possible to make a constraint in EOS independent of 
the gravitational theory. Then, after making a constraint in EOS with the help of observation of the stellar mass, it 
might be possible to restrict on the value of K using the observations of fundamental oscillations, i.e., it might be 
possible to probe the gravitational theory in the strong-held regime observationally. 
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FIG. 4: For stellar models with Madm = IAMq, the frequencies of fundamental and first overtone torsional modes are shown 
as functions of parameter K with EOS A (solid lines with circles) and EOS II (dotted lines with squares). The left panel 
corresponds to the fundamental modes with I = 2, 3, 4, and 5, while the right panel corresponds to the first overtone with 
I = 2. 



V. CONCLUSION 



In this article, focusing on the torsional oscillations inside the crust region, we have derived the perturbation 
equation of neutron stars in TeVeS and calculated their eigenfrequencies. We find that the frequencies of torsional 
modes in TeVeS become larger than those expected in GR, whose dependence on the gravitational theory is stronger 
than that of fluid oscillations involving the variation of density. It is also found that the stellar models with softer 
EOS are more sensitive about the gravitational theory than those with stiffer EOS. Since the frequencies of overtone 
depend strongly on the adopted EOS rather than on the gravitational theory, it is possible, via the observation of the 
torsional frequencies of overtone, to restrict on the EOS independent of the gravitational theory. Then, with the help 
of the observation of stellar mass, it might also be possible to make a constraint in the gravitational theory via the 
observation of fundamental torsional modes. It should be noticed that, since these imprints of TeVeS come from the 
presence of a scalar field, it could be also possible, via observations of such oscillations of neutron star, to probe the 
existence of a scalar field. 
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As a first step, we assume the Cowling approximation in this article, i.e., our examinations are restricted to only 
fluid oscillations, so we should make a more detailed study including the metric and vector field perturbations. Via 
these oscillations, one could obtain additional information, and combining that with the results reported in this article 
would provide more accurate constraints on the gravitational theory in the strong-field regime. Additionally, since it 
is suggested that the nuclear structure in the bottom of a crust could be nonuniform, the so-called "pasta structure," 
we might consider the effect of the pasta phase on the torsional oscillations [23]. Furthermore, we should take into 
account the magnetic effects on the frequencies, because the torsional oscillation could be detected in the magnetars 
as the quasiperiodic oscillations during the giant flares [33| . Considering these additional effects, one must obtain 
more accurate constraints in a more realistic situation. 
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